Noise Predictions for STM in Systems with Local Electron Nematic Order 
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We propose that thermal noise in local stripe orientation should be readily detectable via STM 
on systems in which local stripe orientations are strongly affected by quenched disorder. Stripes, a 
unidirectional, nanoscale modulation of electronic charge, are strongly affected by quenched disorder 
in two-dimensional and quasi-two-dimensional systems. While stripe orientations tend to lock to 
major lattice directions, dopant disorder locally breaks rotational symmetry. In a host crystal with 
otherwise C4 rotational symmetry, stripe orientations in the presence of quenched disorder map to 
the random field Ising model. While the low temperature state of such a system is generally a stripe 
glass in two dimensional or strongly layered systems, as the temperature is raised, stripe orientational 
fluctuations become more prevalent. We propose that these thermally excited fluctuations should 
be readily detectable in scanning tunneling spectroscopy as telegraph noise in the high voltage part 
of the local I(V) curves. We predict the spatial, temporal, and thermal evolution of such noise, 
including the circumstances under which such noise is most likely to be observed. In addition, we 
propose an in-situ test, amenable to any local scanning probe, for assessing whether such noise is 
due to correlated fluctuations rather than independent switchers. 



I. INTRODUCTION 

There is experimental evidence that many strongly cor- 
related electronic systems such as nickelates, cuprates, 
and manganites exhibit some degree of local inhomo- 
gencity, i.e., nanoscale variations in the local electronic 
properties. Describing the electronic behavior of these 
materials involves several degrees of freedom, including 
orbital, spin, charge, and lattice degrees of freedom. Dis- 
order only compounds the problem. Not only can cer- 
tain types of disorder destroy phase transitions, leaving 
mere crossovers in their wake, it can fundamentally al- 
ter ground states, sometimes forbidding long range order. 
Especially in systems where different physical tendencies 
compete, disorder can act as nucleation points for com- 
peting ground states. The interplay between many de- 
grees of freedom, strong correlations, and disorder can 
lead to a hierarchy of length scales over which the re- 
sulting physics must be described. While such electronic 
systems are highly susceptible to pattern formation at 
the nanoscale, unfortunately most of our theoretical and 
experimental tools are designed for understanding and de- 
tecting homogeneous phases of matter. 

Because of these difficulties, the presence of quenched 
disorder in real materials can make it exceedingly dif- 
ficult to discern the character of the local pattern for- 
mation. One might expect to overcome the difficulty by 
growing cleaner samples, but certain forms of local or- 
der (such as stripe orientations, i.e., the nematic com- 
ponent of stripes) are so fragile that any finite disorder 
precludes a long range ordered ground state in two di- 
mensions. (Finite but small disorder can preclude it in 
layered systems.) It is therefore desirable to develop the- 
oretical and experimental techniques aimed at character- 
izing and detecting local electronic order. In previous 
work, we have shown that noise and noncquilibrium ef- 



fects can be used to help illuminate the local electronic 
tendencies^ In this paper, we discuss the implications 
of such noise for scanning tunneling microscopy. In par- 
ticular, we propose that thermally excited switching of 
local stripe orientations should be readily detectable via 
scanning noise microscopy, appearing as telegraph noise 
in the high voltage part of atomic scale dl/dV curves. 

Long-range stripe order has been detected in a sub- 
set of cuprates, where it can coexist with supercon- 
ductivity, as well as in the related nickelates^ How- 
ever, most cuprate superconductors lack evidence of 
long-range stripe order. This calls into question the- 
ories of cuprate superconductivity based on quasi-one- 
dimensional electronic structure. On the other hand, be- 
cause the pairing scale is large, it can in principle be 
established on short length scales. (In fact, STM ex- 
periments have shown that the correlation length of the 
superconducting order parameter, i.e., the distance over 
which the gap scale can vary, is about 10-20A.) Thus the 
salient question is not whether long-range stripe order 
is ubiquitous in cuprate superconductors, but whether 
short-range stripe order is ubiquitous. Nevertheless, the 
fact that long-range ordered stripes have been observed 
in several materials with strong electronic correlations 
means that the study of these structures in the presence 
of quenched disorder is interesting in its own right, re- 
gardless of issues of the mechanism of cuprate supercon- 
ductivity. 

Indeed, recent STM experiments show evidence of an 
electronic cluster glass of locally unidirectional domains 
in Na-CCOC and Dy-Bi2212 at low temperature^ We 
predict that upon raising the temperature in such sys- 
tems, more and more stripe domains will become ther- 
mally excited, resulting in local fluctuations in the stripe 
orientation. The main purpose of this paper is to pre- 
dict the spatial, temporal, and thermal evolution of such 
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noise. 

Telegraph noise in the mesoscopic transport proper- 
ties of a YBCO nanowire has been reported by Bonetti 
et al. in the pseudogap regime^, appearing on timescales 
of about 10-50 seconds, and producing resistance fluctu- 
ations of the order of 0.5% of the total resistance of the 
wire. Wc have previously modeled this telegraph noise 
via a mapping of stripes in the presence of quenched dis- 
order to the random field Ising modeli (as explained be- 
low), along with a further mapping to a related random 
resistor network. While the nanowire transport results 
can be captured by this model, it was not possible to 
definitively rule out other sources of the telegraph noise, 
such as superconducting fluctuations.— Transport noise 
has also been reported on timescales of 1-1000 seconds 
in the c-axis resistivity of bulk LSCO at T < 0.3K by 
Raiccvic and coworkers.— 

Is the low temperature glass of unidirectional domains 
observed in Na-CCOC and Dy-Bi2212 related to the 
transport noise observed in YBCO and LSCO? If sim- 
ilar noise were to be found with STM, it would begin to 
bridge the gap between the real-space imaging capability 
of that probe and the noise seen in mesoscopic and macro- 
scopic systems. This would allow concrete connections to 
be made between the microscopic and macroscopic be- 
havior, and provide guidance as to how to extract infor- 
mation about local order from macroscopic probes when 
microscopic probes such as STM are not feasible. Fur- 
thermore, since transport noise is a bulk measurement, 
finding evidence of similar noise in STM would establish 
a connection between physics happening in the bulk and 
that occurring on the surface. In addition, the wealth of 
information available from STM can be used to set the 
parameters of our model, so that concrete predictions can 
then be made for other measurements. 
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(a) Two Dimensions (b)Three Dimensions 

FIG. 1: Schematic phase diagram of the random field Ising 
model, (a) In two dimensions, the critical disorder strength 
is zero, and for any finite disorder strength the model is in a 
disordered phase at all temperatures, (b) In three dimensions, 
there is a finite critical disorder strength. 

We now describe the mapping of stripe orientations 
in the presence of quenched disorder to the random 
field Ising model.— Stripes in a host crystal such as the 
cupratcs and nickelates tend to lock to favorable lattice 
directions, often either vertically or horizontally along 
Cu-0 (Ni-O) bond directions. This breaks the symmetry 




(a)Color visualization (b)Stripe visualization 



FIG. 2: (Color online) A snapshot of a layered random-field 
Ising model with R = 2.5J, T = 0.5J on a 32 x 32 x 8 cubic 
lattice (see text for details). Only the top layer is shown, 
(a) Here, the two states of the Ising variables, o~i = ±1, are 
represented by red and blue patches, (b) The same results 
visualized as vertical and horizontal stripes. 

of the host crystal from 4-fold to 2-fold rotational sym- 
metry. With two allowed orientations, the nematic com- 
ponent of the stripe order parameter is of Ising character, 
where the Ising variable or "pseudospin" a = ±1 corre- 
sponds to horizontal or vertical stripe patches. The ten- 
dency of neighboring stripe patches to align is modeled 
as a ferromagnetic interaction between nearest-neighbor 
Ising variables. Dopant atoms act like a random field 
on the electronic stripe nematic, so that in any given re- 
gion, the particular dopant arrangement will favor one 
or the other orientation of the stripes locally, producing 
a disorder-dependent pattern of stripe patches with two 
orientations. Thus stripes in a host crystal with dopant 
disorder can be mapped to the random-field Ising model 
(RFIM)j4 

h = - j o-io-j — jj_ o-icrj - y~] (k + h i )o- l (i) 

<*i>|| i 

where J is the in-planc coupling, J± is the interlaycr 
coupling, h is an orienting field, and hi are the random 
pseudofields at each site. In the summations, (ij)u de- 
notes a summation over sites within a plane, and (ij) ± 
denotes summation over sites in neighboring planes. The 
orienting field may be, for example, h oc B x 2 — B y 2 where 
B is the physical external magnetic fields The random 
fields hi are chosen independently from a Gaussian distri- 
bution with standard deviation R (the disorder strength) . 
The order parameter M = ^ ■ (cr.j) represents the de- 
gree of stripe orientation, i.e., nematic order. The size 
of a stripe patch sets the spacing of the Ising lattice, and 
we may use information from experiments to set bounds 
on this lattice spacing. For example, for the purposes 
of comparing with STM experiments on Na-CCOC and 
Dy-Bi2213^, a reasonable lattice spacing would be 16A. 

The materials of interest are layered materials, corre- 
sponding to a quasi-two-dimensional system. As shown 
schematically in Fig. [T] for any finite disorder strength, 
the two-dimensional RFIM remains disordered at all tem- 
peratures in the thermodynamic limit £ In three dimen- 
sions, the RFIM can order below a critical disorder 
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strength. For the finite layered systems we consider here, 
there is a finite critical disorder strength, although it 
may be quite suppressed for strongly layered systems. 
This implies that materials which are strongly layered 
are less likely to support an ordered Ising nematic phase 
than materials with stronger coupling in the c-direction. 
Although the charge component of stripes is subject to 
Coulomb coupling from plane to plane, because they are 
overall charge neutral this coupling is greatly diminished 
for distances longer than the local modulation wave- 
lenth (about 16 A), and the short-range interactions of 
our model are sufficient to describe the c-axis coupling. 



II. THERMALLY ACTIVATED STRIPE 
ORIENTATION SWITCHING 

In Figure [21 we show an example of the low temper- 
ature behavior in the disordered phase, revealing a low 
temperature pattern with locally unidirectional domains 
of various sizes. Fig. Eta) shows a single thermalized 
configuration in our model, and Fig. [5Jb) shows the cor- 
responding stripe visualization. In order to do the cal- 
culation, we have performed Monte Carlo simulations of 
Eqn. [1] with checkerboard Glauber updates^ 7 - for a lay- 
ered RFIM on a 32 x 32 x 8 cubic lattice. We use periodic 
boundary conditions in the x- and ^-directions along with 
open boundary conditions in the z-direction, and we show 
results from the top layer, as a rough approximation to 
the geometry of an STM experiment on a surface. We 
have used coupling constants appropriate to a layered 
system, J = 1 within each layer and J z = 0.125 be- 
tween layers, with a disorder strength R = 2.5 J, which is 
above the critical disorder strength. This typical config- 
uration is taken after approximately 10 7 timesteps. (It 
takes about 5 x 10 4 steps to thermalize at these param- 



eters.) In Figure 2(a) we show a conventional visualiza- 



tion of the Ising state, using color to represent the state 
of the spin. In Figure 2(b) we have shown the director 



of the Ising nematic, in order to represent the average 
direction in which stripes run at any given site in the 
corresponding Ising model. 

Because we have used a disorder strength which is 
larger than the critical disorder strength, the ground 
state configuration is dominated by quenched disorder, 
and although there are domains of various sizes, the av- 
erage domain size is small. A consequence of this is 
that upon thermal cycling, the low temperature state 
returns to the same pattern of stripe orientations, since 
the local orientation of each stripe patch is dominated 
by the quenched disorder rather than by the thermal his- 
tory. The tunneling asymmetry maps of Ref. bear a 
strong resemblance to configurations of the RFIM (see 
Fig. 2(b) |. Note that in Fig. [5] we have shown a single 
(typical) configuration, rather than an ensemble average. 
We believe that this is the correct comparison for STM 
experiments at T = 4.2K (as opposed to an ensemble 
average), since the timcscalcs over which a single dl /dV 



curve is recorded arc small compared to the switching 
time of stripes at that temperature^ 

While most fluctuations are frozen out at low temper- 
ature, upon raising the temperature, certain regions of 
the system become active, and exhibit thermally driven 
stripe orientation fluctuations. In an STM experiment, 
this should lead to local telegraph noise in the active 
regions, where some atomic positions will alternate be- 
tween having a charge stripe on them or not. We show 
below where to look for such noise as a function of posi- 
tion and temperature, and then describe an algorithm for 
detecting whether such noise is due to correlation effects 
rather than from some local, noninteracting source. 

Note that microscopically, a rigid rotation is not re- 
quired in order to produce the thermally fluctuating lo- 
cal Ising variable of our model. Rather, minor local rear- 
rangements are sufficient to change the direction of any 
particular unidirectional domain. Exactly how such a lo- 
cal rearrangement transpires is a question of the local en- 
ergy barriers of a disordered, many-body system, and will 
be dominated by whatever is the lowest local activation 
energy in any given region. Our proposal is independent 
of these microscopic details. It is, of course, conceivable 
that other mechanisms (e.g., side-to-side stripe fluctua- 
tions in a disordered system or fluctuations of other in- 
homogeneous patterns) may give rise to fluctuating Ising 
variables. 

Because the materials in question are extremely com- 
plicated and each Ising variable in our model encom- 
passes many electrons, it is difficult to arrive at a re- 
liable microscopic estimate of timescales. Nevertheless, 
here is a rough argument to show that the numbers are 
not unreasonable. Assuming that the stripe fluctuations 
involve electron dynamics, the attempt frequency would 
be of the order of the Fermi velocity (2 eVA) divided by 
the stripe spacing (16 A), that is, 10 15 Hz. Noise in the 
YBCO nanowire was observed on a timescale of about 50 
seconds at 100K4 For an Arrhenius type law, this cor- 
responds to a local energy barrier of about 3800K, not 
far from the bare Hciscnberg spin exchange couplings in 
the cuprates. While the Heisenberg coupling is not equal 
to our Ising coupling, it is one parameter which sets the 
Ising coupling, and so the two may be of the same mag- 
nitude. Note, however, that due to the randomness in 
our model, every cluster has a different timescale, many 
clusters have multiple timescales, and these timescales 
span several orders of magnitude. As a result, the issue 
is less one of whether timescales will be accessible exper- 
imentally, and more one of where spatially to look for 
fluctuating domains. We address this issue in Fig. [3] 

Figure [3] shows what happens when the temperature 
is raised on the same system as that in Figure [5] In the 
leftmost column of Fig. [31 the grey and black regions are 
exhibiting stripe orientation fluctuations. In the corre- 
sponding stripe visualization (second column of Fig. [3]) , 
this appears as places where the different stripe orienta- 
tions appear to "invade" each other. These active regions 
have relatively high variance compared to the surround- 
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FIG. 3: (Color online) Fluctuations of the RFIM for a disorder strength of R = 2.5J at various temperatures as noted in 
the figure, on a 32 x 32 x 8 cubic lattice (see text for details). Only the top layer is shown. In the leftmost column, we show 
the time-averaged (pseudo)spins (en), where blue and red denote positive and negative averages, respectively, with shades of 
grey denoting average pseudospins near zero. In the next column, we plot each averaged pseudospin variable as the director 
of the Ising nematic, to visualize the pattern of horizontal and vertical stripes. Red sites from the first column correspond to 
horizontal stripe orientation, and blue sites from the first column correspond to vertical stripe orientation. The third column 
shows the variance of the (pseudo)spins, ((A<7i) 2 ). The rightmost column shows the local switching timescale n (see text for 
definition). At low temperatures, fluctuations are all but frozen out, and local switching timescales diverge. As the temperature 
is raised, more and more clusters become active, showing as bright spots in the variance maps (third column). Note that the 
boundaries of the low temperature clusters become active first as temperature is raised. At high temperature, noise is apparent 
throughout the sample. 



ing regions, as shown in the third column of Fig. [3] The 
black patches in the leftmost column can be seen as active 
(white) patches in the variance maps. As temperature is 
raised, more and more regions in the sample become "ac- 
tive". The places which become active first are bound- 
aries between the low temperature (frozen) domains, and 
especially domain boundaries with several small domains. 
These are the best places to look for local telegraph noise 
upon raising temperature. 



Of course, the observation of telegraph noise is not 
sufficient to determine whether correlations are present 
in the system. Telegraph noise could arise from other 
sources, such as an unstable defect on the tip, or even 
defect hopping on the surface of the material itself. How- 
ever, these are uncorrelated sources of noise due to iso- 
lated, independent switchers. We are interested in how 
to identify noise arising from a correlated system. One 
way to distinguish such "spurious" possible sources of 
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FIG. 4: (Color online) A line of sites in the RFIM at R = 2.5 J 
and T = 0.5J. This corresponds to the simulation from the 
fourth row of Fig. [3] and we have chosen a line of sites going 
through an "active" region. The colored arrows in the inset 
correspond to the colors of the time traces shown in Fig. [5] 
and of the block variances shown in Fig. [8] 



telegraph noise from that due to correlations in an inter- 
acting model such as the one considered here is to ana- 
lyze the power spectrum of a time series of the data. It 
is well known that for a "stationary" random telegraph 
signal (i.e., one with a constant probability of flipping 
from "up" to "down" , and with a separate unchanging 
probability of flipping from "down" to "up") the result- 
ing power spectrum is a Lorentzian;^ and deviations 
from this form would presumably indicate correlations^ 
However in the materials of interest, where the timcscalcs 
have been observed to be as long as 10-50 seconds in 
YBCO^ and 1000 seconds in LBCO^, a time domain cri- 
terion is much more desirable. 

We first describe what we mean by correlated switch- 
ers, and then we present a criterion for detecting spa- 
tial correlations using a time domain analysis at a single 
site. In Figure HJ we show a few sites in a larger system, 
with color-coded arrows denoting particular local stripe 
patches. We have chosen some sites which are mostly 
"vertical" (blue), and some which are mostly "horizon- 
tal" (red), along with some sites which are active switch- 
ers (grey). The colors of the arrows denoting particular 
sites in Fig. 2] correspond to the colors of the time traces 
of Ui in Fig. [5j The local orientation at each site displays 
telegraph noise, but the nature of the telegraph noise is 
different for different sites. For example, the red trace 
shows a few sharp spikes, whereas the blue and green 
traces show frequent switching. Notice that the noise is 
correlated: when the blue spin is "up" , the green spin is 
also up; however, when the blue spin is down, the green 
spin becomes "active" and switches rapidly. This behav- 
ior indicates correlated switchers, which is the essence 
of the interacting model we use. One might imagine 
that such correlations could be detected via simultaneous 
measurements with two (or more) STM tips; however, 
that would require placing the tips closer together than 
can be achieved with current setups. We therefore focus 
on developing a criterion for detecting correlated noise 
with a single tip. 

One way to characterize the temporal behavior of 
the local noise is to study the autocorrelation function, 
^i{t) ~ (ci(0)oj(4)}. However, we have found that the 
autocorrelation is prone to large statistical errors at long 



times, and hence provides an unreliable characterization 
of the data. The origin of this issue lies in the interac- 
tions of the system. Because each site is interacting with 
several other sites, each site responds at several different 
timescales. The true asymptotic behavior of the local 
autocorrelation is not achieved until well beyond all of 
these timescales. 

For this reason, we instead use the "on-the-fly" re- 
blocking algorithm 9 - described below, which is efficient, 
robust, and lends itself easily to real-time signal process- 
ing during data collection. A graphical representation of 
this method is shown in Fig. [5] A block spin is defined to 
be the average spin at a particular site, averaged over m 
consecutive timesteps, where m is the block size. As the 
Monte Carlo simulation proceeds, we update the values 
of the most recent complete blocks of sizes that are pow- 
ers of two: m = 1,2, 4, 8, 16, ... . We also accumulate the 
sums and sums-of-squares of block spins, so that we can 
ultimately construct the block averages and block vari- 
ances. We use an unbiased estimator for the variance, 
so that for, e.g., a 16-element time series at block size 
m — 2, 
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This procedure requires minimal computational effort, 
and furthermore requires the storage of only C(log 2 A r ) 
numbers in memory, making the method amenable to 
real-time data analysis. The variance ^(A<r) 2 ) m as a 
function of the temporal block size m can then be used 
to distinguish a telegraph signal that is free of spatial 
correlations (i.e., one due to a stationary random tele- 
graph signal) from a telegraph signal due to an interact- 
ing model, the spatial correlations of which contribute 
multiple timescales to each local switcher. 

We first describe the results of this analysis for an un- 
correlated, independent switcher (a stationary random 
telegraph signal (RTS)), and then contrast it with local 
telegraph noise from an interacting model, which encodes 
information about the correlations. At the block size m, 
the block variance can be expressed as 
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where the average (• • • ) is taken over all blocks of size m, 
and ip(t) is the time-dependent part of the autocorrela- 
tion function Vt'(t). For a stationary discrete-time RTS§., 
the autocorrelation is exponential (corresponding to a 
Lorentzian power spectrum): 



where 
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1/T = X = l-p-q, (6) 
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FIG. 5: (Color online) Time traces of the nine sites denoted in Fig. [4] with R = 2.5J and T = 0.5J. The colors of the time 
traces correspond to the colors of the arrows in Fig. [4] 
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FIG. 6: (Color online) Visualization of the "real-time" reblocking algorithm for the spin history at a single site. The horizontal 
axis represents Monte Carlo time (or real time in a discrete time series from experiment); the vertical axis represents the block 
size, m. Dark areas represent previously computed blockspin values that have already have been discarded from memory. 
White areas are yet to be calculated. The lightly shaded (green) areas denote the current state which must be held in memory, 
of order 01og2./V at time step t = N. 
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FIG. 7: (Color online) Normalized block variances 
((Acr( m )) 2 ) / ((Act^)) 2 ) vs. block size m of an uncorrelated 
switcher (stationary RTS) for various characteristic timescales 
r, showing the crossover from constant to 1/m behavior. 



where p (q) is the conditional probability of switching 
from the "up" ( "down" ) state to the "down" ( "up" ) state 
at the next time step, and r is the resulting characteristic 
timescale. Thus the block variance is 



((Aa (m) ) 2 ) =a 



m - m\ 2 + 2A(A m - 1) 
m 2 (A-l) 2 ' 



(7) 



Figure [7] shows the normalized block variance 
((A<T( rn )) 2 ) /((AtT(i)) 2 ) vs. m for uncorrelated switchers 
with various r. The small m behavior is a horizontal 
line on a log- log plot, whereas the large m behavior is a 
straight line corresponding to ((Acr( m )) 2 ) oc 1/m. The 
extrapolated intersection of these two linear regimes 
may be used to extract a rough timescale r 2 ; for each site 
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where m max is the largest block size for which variance 
data is available. We will see that deviations from the 
horizontal line at low values of m can be used as a local 
measure of spatial correlations. 

In Fig. [SI we show the variance in our model, where 
spatial correlations have caused a deviation from the sim- 
ple behavior of a stationary RTS. The colors of the curves 
correspond to the colors of the arrows in Fig. |4j A couple 
of these curves mimic that of a stationary RTS. For ex- 
ample, the red curve begins with the horizontal line char- 
acteristic of a stationary RTS, and then crosses over to 
the 1/m behavior at a rather large timescale. In addition, 
the purple curve displays the 1/m behavior throughout 
the plot, consistent with a local stationary RTS with a 
rather short timescale. The rest of the curves in Fig. [5] 
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FIG. 8: (Color online) Block variances V t vs. block size 
m for the line of sites shown in Fig. [4] The colors of the 
block variance curves correspond to the colors of the arrows 
in Fig. [4] For the red and green traces, the intersection of 
the dotted lines shows the characteristic timescale of those 
sites, Ti, as described in the text. While the results for the 
red and purple traces mimic that of a stationary RTS with 
a long and short timescale, respectively, the other timetraces 
show strong deviation from this behavior, a clear signature of 
the correlations in the system. There are no brown and grey 
traces because the corresponding spins have zero variance over 
the sampled time interval. 



show a striking deviation from the stationary RTS be- 
havior, in that they begin at small m with a slope which 
is in between a horizontal line and 1/m behavior. This 
may be taken as an indication of correlations. In sev- 
eral of the curves, multiple "knees" are evident within 
the same curve, further evidence of correlations, in that 
there are multiple timescales present in the behavior of a 
single site. This type of analysis is amenable to real-time 
data processing of a telegraph STM signal in an "active" 
region, in order to determine whether the noise is due 
to a local independent switcher (which may be due to, 
e.g., a "dirt" effect), or whether it is due to the more in- 
teresting case of correlations due to interactions among 
the electronic degrees of freedom^ In addition, it is free 
from the windowing errors that an FFT would entail. 

In the right column of Fig. [3j we plot maps of the lo- 
cal switching timescale, r, estimated from Eqn.[5] These 
may be viewed as a map of the dynamic clusters, i.e., 
those sites which fluctuate together. Note that this is 
different from the identification of the low temperature 
domains. The behavior of the largest timescales as tem- 
perature is lowered (i.e., in the active regions) is consis- 
tent with an Arrhenius law. However, a single timescale 
is insufficient to characterize the local dynamics because 
of the correlations, and there arc longer timescales in the 
active regions than can be properly represented with a 
2D color plot. 



III. DISCUSSION 

We have argued that there should be telegraph noise 
in STM experiments on locally striped materials at in- 
termediate temperatures. Our proposal can be summa- 
rized as follows: Using STM in, e.g., a tunneling asymme- 
try (TA) mode, the low temperature spatial map can be 
used to identify the low temperature domain structure. 
Then, place the tip on a domain boundary, preferably 
where several small domains arc clustered, and take time 
traces of the signal at a few sites. If no telegraph noise 
is evident, raise the temperature, and repeat. At high 
enough temperature, switching of the stripe orientation 
should become thermally excited, resulting in local tele- 
graph noise. We have furthermore proposed a criterion 
by which such noise can be analyzed in the time domain, 
in order to determine whether it is arising from a local, in- 
dependent switcher, or whether it is instead arising from 
correlated fluctuations in an interacting model. 

Important issues include which signal should be 
tracked, at what timescale do we expect telegraph noise 
to be evident, and what are good candidate materi- 
als for testing our proposal. Some evidence of stripe 
structure has been reported in STM on BSCCQi ^ 
in the 8-15meV range, derived from measurements of 
dl/dV. However, spatial modulations in the low en- 
ergy response have also been attributed to quasiparticlc 
interference! 1 ^ 13 More dramatic evidence of locally uni- 
directional electronic domains, with minimal signal pro- 
cessing involved, has been reported in Dy-BSCCO and 
NCCOC2 via tunneling asymmetry (TA) maps. Specif- 
ically, the ratio of the tunneling current at ±150mV re- 
sulted in high contrast images of domains with local Ising 
symmetry, consistent with the model studied here. For 
the purposes of studying noise associated with these Ising 
domains, a time trace of the current near ±150mV would 
be a good place to start. (Note that the magnitude of 
variation in signal in those experiments is on the order of 
30%.) The variance test proposed here, which depends 
on the time domain properties of a discrete signal change, 
would apply whether the "tip height" z was held constant 
and a time trace of the current / measured, or whether 
the current I was held fixed and a time trace of z was 
recorded. 

It is worth emphasizing that telegraph noise has al- 
ready been reported in the transport properties of a 
YBCO nanowirei 3 - In addition, similar behavior was 
noted in a BSCCO nanowir&i^, although sufficient data 
was not gathered to demonstrate telegraph noise in a 
time trace. In this case, switches in the resistivity were 
observed starting above T c (since the resistivity is zero 
below T c ), and going as high as T = 150/\ . The noise was 
smaller in magnitude, and less frequent than the corre- 
sponding noise in YBCO. However, BSCCO is much more 
amenable to an STM experiment, and given that locally 
unidirectional domains have been observed in TA maps 
in Dy-BSCCO, this system is a good candidate for ob- 
serving thermally excited stripe switching noise via STM. 
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In this paper we argue that fluctuations of stripe ori- 
entations will produce correlated telegraph noise in scan- 
ning tunneling experiments, and we propose a method 
for analyzing such noise in the time domain. This does 
not, of course, rule out the possibility of correlated noise 
emerging from other sources (e.g., side-to-side stripe fluc- 
tuations, or fluctuations of other inhomogeneous pat- 
terns, such as checkerboards or quasi-particlc interference 
patterns). 

IV. CONCLUSIONS 

In conclusion, we have analyzed the local noise char- 
acteristics of the random field Ising model in the dis- 
ordered phase. For materials with local stripe corre- 
lations in the presence of disorder, this model predicts 
that there will be thermally excited fluctuations in the 
stripe orientation. We propose that such fluctuations 
should result in correlated telegraph noise in the STM 
response near stripe domain boundaries as temperature 



is raised. The use of this method in conjunction with 
TA scans can establish whether there is a connection be- 
tween the observed low-T unidirectional, glassy domains 
in Na-CCOC and Dy-Bi2212, and the transport noise 
observed in YBCO and LSCO, and also put constraints 
on our model. We have furthermore proposed that the 
block variance of telegraph noise can be used as a purely 
local indicator of spatial correlations. The criterion can 
be applied to any local probe (AFM, MFM, etc.) and 
because it is accumulated directly in the time domain, it 
is amenable to real-time signal processing. 
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Other examples of an external perturbation which breaks 
the rotational symmetry of the lattice and would there- 
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In some magnetic systems Kawasaki dynamics may be 
more appropriate, as they conserve local magnetization. 
However, in our problem there is no reason for nematic- 
ity to be locally conserved, so Glauber dynamics are ap- 
propriate. Furthermore, because the switching of the local 
nematic orientation (a single Ising variable in our coarse- 
grained model) involves many electrons, we do not expect 
this to be a coherent quantum process. Therefore, we treat 
the dynamics classically, rather than using a Caldeira- 
Leggett model. 

Recall that at 100K, telegraph noise in the YBCO nanowire 
has a timescale of 10-50 seconds. By 4.2K, thermal fluc- 
tuations in stripe orientation should be frozen out in that 
material. Likewise, fluctuations in stripe orientation should 
be similarly frozen out in BSCCO and NCCOC at 4.2K 
Second spectra, another frequency-domain criterion, may 
also be used.— 

Note that we have not focused on critical behavior or on 
scaling behavior.— ^ Comparing those type of predictions 
to STM requires that the measurements be averaged over a 
large field of view. One advantage of our proposal is the fo- 
cus on local physics, without the need for spatial averages. 



